Quantum Tunneling of Normal-Superconducting Interfaces in a Type-I 

Superconductor 
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Evidence of a non-thermal magnetic relaxation in the intermediate state of a type-I superconducor 
is presented. It is attributed to quantum tunneling of interfaces separating normal and supercon- 
ducting regions. Tunneling barriers are estimated and temperature of the crossover from thermal 
to quantum regime is obtained from Caldeira-Leggett theory. Comparison between theory and ex- 
periment points to tunneling of interface segments of size comparable to the coherence length, by 
steps of order one nanometer. 

PACS numbers: 74.25. Ha, 74.50. +r, 75.45.+J 

Quantum tunneling of relatively macroscopic solid- 
state objects like flux lines in type-II superconductors 
[l|, [2| and domain walls in magnets [3| have been sub- 
ject of intensive research in the past. Pinning by defects 
and impurities creates complex potential landscape that 
traps flux lines and domain walls inside metastable en- 
ergy minima. Their interaction with environment makes 
this problem the one of macroscopic quantum tunneling 
with dissipation [4J . The latter is especially important for 
the tunneling of flux lines because of their predominantly 
dissipative dynamics [5|-|9( . If individual pinning centers 
are small compared to the dimensions of the vortex or 
the width of the domain wall, the pinning is collective. 
In this case the energy barriers and spatial scales of ther- 
mal and quantum diffusion of flux lines and domain walls 
are non-trivially determined b y s tatistical mechanics of 
the pinning potential [l|, Q, LL1A LLL [ ■ 




FIG. 1: Color online: Schematic view of the interface be- 
tween normal and superconducting regions of type-I super- 
conductor, slightly perturbed by randomly scattered pinning 
centers. 



When placed in the magnetic field, type-I supercon- 
ductors do not develop flux lines. Instead, they exhibit 
intermediate state in which the sample splits into normal 
and superconducting regions separated by planar inter- 
faces of positive energy [12h14 | . Recently, there has been 



a renewed interest to the equilibrium structure, pinning, 
and dynamics of interfaces in type-I superconductors |15l - 
[20( . Pure lead has been mostly used as a prototypical ex- 
perimental system. In the presence of pinning centers the 
interfaces adjust to the pinning potential by developing 
curvature as is schematically shown in Fig. [T] Pinning 
by point or small- volume defects should result in a broad 
distribution of energy barriers. It is, therefore, plausi- 
ble that at low temperature type-I superconductors con- 
tinue to relax towards equilibrium via quantum diffusion 
of interfaces. This situation is similar to the diffusion of 
domain walls in disordered ferromagnets with one essen- 
tial difference. Contrary to a ferromagnetic domain wall, 
the dynamics of the planar interface in a superconductor 



should be dominated by dissipation. 

At low temperature the decay of metastable states cre- 
ated by pinning provides slow relaxation of magnets and 
superconductors towards thermal equilibrium. This re- 
laxation is known as magnetic after-effect. At finite tem- 
perature it may occur via thermal activation with a prob- 
ability proportional to exp(— Ub/T) where Ub is the en- 
ergy barrier. As T — > thermal processes die out and 
the only channel of escape from the metastable state be- 
comes underbarrier quantum tunneling. Its probability is 
proportional to exp(—I e ff/K) where J e // is the effective 
action associated with tunneling. The pre-exponential 
factors in the two expressions are of lesser importance 
because the dependence of the probability on the param- 
eters is dominated by the exponents. Equating the two 
exponents, one finds that the crossover from thermal ac- 
tivation to quantum tunneling occurs at Tq sa hUs/Ief}- 
Experimental evidence of such a crossover in type-II su- 
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FIG. 2: Color online: Magnetic relaxation at various temper- 
atures. Logarithmic time dependence provides an accurate fit 
to the data. 



FIG. 3: Temperature dependence of the magnetic viscosity. 
S(T) tends to a non-zero temperature-independent value in 
the limit of T -+ 0. 



perconductors has been overwhelming [21]. There has 
also been sonic experimental evidence of quantum dif- 
fusion of domain walls in disordered ferromagnets 22 1. 



However, to our knowledge, no literature exists on non- 
thermal magnetic relaxation in type-I superconducors. 
Experimental evidence of such a relaxation and its theo- 
retical treatment are subjects of this Letter. 

The lead sample used in our experiments was an octag- 
onal disk of thickness 0.2 mm and surface area 40 mm 2 , 
cut from a commercial Pb rod of purity 99.999%. It was 
annealed during one hour at 290°C in glycerol and nitro- 
gen atmosphere to reduce the mechanical stress from de- 
fects that might have been introduced during preparation 
of the sample. Magnetic measurements were performed 
with the use of a commercial superconducting quantum 
interference device (SQUID) magnetometer in the field 
up to 1 kOc in the temperature range 2K - 8K with ther- 
mal stability better than 0.01 K. Isothermal magneti- 
zation curves were measured to obtain the temperature 
dependence of the thermodynamic critical field. The fit 
of the data by B C (T) = B C (0)[1 - (T/T c ) 2 } produced 
B c (0) = 802 ± 2 0e and T c = 7.23 ± 0.02 K, in accor- 
dance with the values of the critical field and transition 
temperature reported for lead. 

Magnetic relaxation was measured by first applying the 
field B > B C (T), then subsequently switching the field off 
and recording isothermal temporal evolution of the rem- 
nant magnetization M iem (T) in a zero field. Fig.[2]shows 
the time evolution of M rem (t)/M rem (0) between 2.00 K 
and 7.00 K in steps of 0.50 K. At all temperatures the 
observed slow relaxation followed very well the logarith- 
mic time dependence, M rcm (t) = M rcm (0)[l — S(T)lnt], 
where S(T) is the so-called magnetic viscosity. Tempera- 
ture dependence of S(T) is shown in Fig. [3J Remarkably 
it does not extrapolate to zero in the limit of T — > but, 



instead, tends to a finite temperature-independent limit 
as the sample is cooled down. 

As is well known, logarithmic time-dependence of the 
relaxation is an indication of the broad distribution of 
the energy barriers for the escape from metastable states 
3]. The finite value of S(0) points towards quantum 
mechanism of the escape. By analogy with type-II su- 
perconductors, where non-thermal magnetic relaxation 
is due to quantum tunneling of flux lines, it is reason- 
able to assume that in type-I superconductors the effect 
is due to quantum tunneling of interfaces separating nor- 
mal and superconducting regions. The structure of the 
interface (see Fig. 0]) is determined by two parameters: 
the coherence length £ and the London length Al . Type- 
I superconductivity corresponds to k = A^/C < l/v2- 
Concentration of Cooper pairs l^l 2 gradually goes to zero 
on a distance £ as one moves through the interface from 
the superconducting to the normal region. When cross- 
ing the interface in the opposite direction one would sec 
the magnetic field going down from its thermodynamic 
critical value B c to zero on a distance S = \/\l,£ < £ ■ 

The energy of the unit area of the interface is [23] 
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(1) 



Pinning provides curvature of the interface, see Fig. [T] 
We shall describe such an interface by a singled-valued 
function Z(x, y, r). The potential energy of the interface 
consists of two parts: 
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FIG. 4: Structure of the interface between normal and super- 
conducting regions of type-I superconductor. The magnetic 
field decays on a scale 8 — \/\l£,, while the modulus of the 
Cooper-pair condensate wave function changes on a scale £. 



The first integral is the energy of a two-dimensional elas- 
tic manifold. The second integral is the energy due to 
the pinning potential U(x,y,z). Metastable equilibrium 
is achieved through balance of the elastic energy and pin- 
ning energy that corresponds to the minimum of Eq. (J2J . 
We shall assume that magnetic relaxation occurs due to 
the decay (or formation) of the bumps in the interface 
shown in Fig. [T] We shall describe such a bump by the 
lateral size L and height a. For a particular bump these 
parameters are determined by the local pinning potential. 
Since the latter is unknown we shall test self-consistency 
of the approach based upon theory of tunneling with dis- 
sipation [4j] by extracting the average values of L and a 
from experiment. 

Let us first estimate the energy barrier associated with 
the bump. It is easy to see that the change in the elas- 
tic energy of the interface due to formation of the bump 
(see Fig. [5]) is independent of L and is generally of or- 
der ana 2 . (This follows from the fact that the area of 
a spherical segment above any cross-section of a sphere 
differs from the area of that cross-section by ira 2 ). This 
energy must be balanced by the negative energy of the 
pinning to make the bump an equilibrium state of the 
interface. Consequently, 
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(3) 



with the average value of a should represent the typical 
amplitude of the random pinning potential. 

We want to find the WKB exponent, I e ff/h, for the 
tunneling of Z{x, y) between two configurations of the 
interface corresponding to the local energy minima (see 
Fig. [5]). Same as for the flux lines [5H9J we shall assume 
that the tunneling probability is dominated by the dissi- 
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FIG. 5: Flattening (or formation) of a bump via quantum 
tunneling of a pinned interface (I) separating normal (N) and 
superconducting (SC) regions. 



pation part of the Caldeira-Leggett effective action [4j : 
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(4) 



where r\ is a viscous drag coefficient describing dissipative 
motion of the interface and r = it is imaginary time. For 
a segment of the interface of size L, that tunnels by a 
distance a, the T — value of the effective action in Eq. 
(|4]) can be estimated as 
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(5) 



The drag coefficient r\ can be obtained from the argu- 
ment similar to that of Bardeen and Stephen for the flux 
lines [24j . Let magnetic field be in the y-direction. In the 
presence of the current of density j in the x-direction, the 
magnetic force experienced by the dxdy element of the 
interface in the z-direction is 



dF = — j dxdydz jB . 



(6) 



Writing j via the electric field and normal-state resistiv- 
ity p n as j = E/p n , and substituting here E = (V/c)B 
for the electric field produced inside the interface moving 
at a speed V in the z-direction, one has j — (V/c)(B/p n ). 
This gives 



dF 
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for the force per unit area of the interface. Substitution 
into this formula of B w B c exp(— z/5) finally yields 



dF 
dxdy 



T)V, 



V 



2p v 



(8) 



As has been explained in the introduction, the 
crossover from thermal to quantum diffusion of the in- 
terface should occur around Tq = hUB/Ieff- With the 
help of Eqs. ©, ©, and © one obtains 



Tr, 



4ir 2 ~h<7 4ir\/2hp rl 
rjL 2 3^/kL 2 



(9) 



Notice that due to the dimensionality of the problem 
Tq does not depend on the size of the tunneling step 
a. Recalling the expression for Xl in terms of the ef- 
fective mass m and concentration n of the electrons, 
Xl = [?7ic 2 /(47re 2 n)] 1 / 2 , and writing p n = (mv/e 2 n) — 
4irisX L /c 2 in terms of the normal electron collision fre- 
quency v, the crossover temperature can be presented in 
the form 
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that shows its explicit dependence on the microscopic 
parameters of the material. 

Let us now compare our experimental findings with 
theoretical results. The temperature of the crossover 
from thermal activation to quantum tunneling can be es- 
timated from the following argument. At non-zero tem- 
perature the magnetic viscosity S shown in Fig. 3 has 
contributions from both, thermal activation and quan- 
tum tunneling, S = St + Sq, where Sq — S(0). The 
parameter Tq is defined as temperature at which the 
two contributions are equal, that is, St — Sq and 
S(T Q ) = 2S Q . This gives T Q in the ballpark of 4K. 

The values of Xl and £ in lead are 37nm and 83nm, 
respectively, giving k = Xl/S, = 0.45. Thermodynamic 
critical field, B c , is close to 800G. For the energy of the 
unit area of the interface Eq. ([T]) gives a ~ 0.4 erg/cm 2 . 
Normal resistivity of lead at 4 K is of order 25[ 5 x 10~ u 



Q- m w 5.6 x 10 



-21, 



Eq. 



(|8J) then gives for the drag 
coefficient r/ ca 0.35 erg- s/cm 4 . We shall now check self- 
consistency of our model by computing the average size 
of the tunneling segment L and the tunneling step a. 
From Eq. ([9]) one obtains L w 90 nm ~ £, which is rather 
plausible. Indeed, L ~ £ describes the segment of the 
interface inside which Cooper pairs are strongly corre- 
lated and, therefore, they can collectively participate in 
a coherent tunneling event. For the tunneling transition 
to occur in our experimental time window of one hour, 
I e ff cannot significantly exceed 2bh. According to Eq. 
([5]) this condition is satisfied by tunneling steps a below 
1 nm, which is also quite plausible. According to Eq. ([3]). 
the typical energy barrier must be then of order 100K in 
accordance with the fact that thermal activation dies out 
below 4K. 

In Conclusion, we have observed non-thermal magnetic 
relaxation in lead that we attribute to quantum tunnel- 
ing of small segments of interfaces separating normal and 
superconducting regions. Theory of such a tunneling has 
been developed. Comparison between theory and exper- 
iment suggests macroscopic quantum tunneling of inter- 
face segments comparable in size to the coherence length, 
by steps of order one nanometer. 
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